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ABSTRACT 

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the 
Crab nebula lasting for approximately nine days in April of 20f 1. The source, which at optical 
wavelengths has a size of «11 ly across, doubled its gamma-ray flux within eight hours. The 
peak photon flux was (186 ±6) x 10 -7 cm -2 s _1 above 100 MeV, which corresponds to a 30-fold 
increase compared to the average value. During the flare, a new component emerged in the 
spectral energy distribution, which peaked at an energy of (375 ± 26) MeV at flare maximum. 
The observations imply that the emission region was likely relativistically beamed toward us and 
that variations in its motion are responsible for the observed spectral variability. 
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1. Introduction 

The Crab nebula is the remnant of a supernova 
observed in 1054 AD. The explosion left behind 
a rotating neutron star emitting electromagnetic 
radiation pulsed at the rotation period, that pow- 
ers a wind of relativistic particles. These particles 
interact with the remnant gas and magnetic field, 
causing the nebula to glow brightly at all wave- 
lengths, predominantly by synchrotron radiation. 
The spectral energy distribution (SED) of the neb- 
ula is, accordingly, dominated by a synchrotron 
component extending from radio wav elengths into 
the gamma-ray band ( Hester 2008f) . Above 450 



MeV a second component emerges, attributed to 
inverse- Comp ton scattering by the same relativis- 



tic particles (IGould Il965c Ide Jager fc Harding . 
19921 lAtovan et al. lll996l) . The angular size of 
the Crab nebula is «0.1° in the optical and smaller 
at higher energies. This corresponds to 3.5 pc, or 
11 ly, at its estimated distance of 2 kpc ([Trimble 
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19731) . 

Today, the pulsar and nebula (henceforth re- 
ferred to together as the Crab) are considered 
prime examples of non-thermal sources in the Uni- 
verse and serve as a laboratory for relativistic 
plasma physics. New puzzles for our understand- 
ing of the Crab have been posed by the detection 
of three bright gamma-ray flares by the AGILE 
satellite and the Large Area Telescope (LAT) on 
board the Fermi sat e llite between 2007 a nd 2010 
(|Abdo et al. I l2011at iTavani et al~l l201l[) . Dur- 
ing these flares the unpulsed component of the 
gamma-ray flux increased by a fact or of »10 o n 
time scales as short as 12 hours (Balboetal 



20111 ). while the period and flux of the pulsed com- 
ponent remained stable. 

More recently, in April of 2011 the LAT de- 
tected a fourth flare, three times brighter than 



any of the previo us ones (jBuehler et al. 1 12011 
Havs et al. 1 12011 ). The flare wa s swiftly con- 
firmed by the AGILE satellite (|Striani et al. 



20111 ). Observations at lower energies, in partic- 
ular by the Chandra X-ray observatory, have not 
yet revealed any variability correlated with the 



gamma-ray flares (jTennant et al. 1 120111 ) . How- 



ever, analysis of these observations is ongoing and 
will be discussed elsewhere. Here we present the 
LAT gamma-ray results obtained during the flare, 
together with a broader analysis of the first 35 
months of Crab observations by Fermi. 

2. The Large Area Telescope and data 
analysis 

The LAT is a pair-conversion telescope, sensi- 
tive to gamma rays with energies greater than 20 
MeV. It has a large field of view («2.4 sr) and 
images the full sky every three hours. The angu- 
lar resolution of the LAT varies with photon en- 
ergy. The 68% containment radius ranges from 
appro ximately 6° at 70 M eV to 0.2° above 10 
GeV (|Atwood et al. 1120091) . LAT scientific obser- 



vations began in August 2008. We analyzed data 
taken within 20° of the Crab in the first 35 months 
of observations (MJD 54683-55728). The average 
spectral properties of the Crab were derived from 
the first 33 months of observations (MJD 54683- 
55664), excluding the April 2011 flare. 

Fluxes and spectra were obtained by maxi- 
mizing the likelihood of source models using un- 



binned gtlike from the Fermi Science Tools 9- 
23-01. The models included all sources in the 
second LAT s ource catalog withi n 20° of the 



Crab position (jAbdo et al. 1 1201 lq) plus models 



for the Galactic and isotropic diffuse emission 
(gal_2yearp7v6_v0 , iso_p7v6source). The pa- 
rameters left free to vary in the likelihood fit 
for 33-month average spectra were the spectral 
parameters of the Crab, the normalization of 
the diffuse components and a power-law spec- 
tral index for scaling the Galactic diffuse emis- 
sion model. For analysis on shorter time scales 
only the isotropic diffuse normalization was var- 
ied along with the Crab spectral properties. All 
other parameters were fixed to the 33-month av- 
erage maximum likelihood values. The Sun was 
included in the source model when it was within 
20° of the Crab. The solar spectra during the 
two passages in fron t of the Crab were taken from 



Abdo et al. I (|2011bl ). The Moon was not included 
in the source model as its gamma-ray contamina- 
tion was found to be negligible for observations 
presented here. 

We used the P7_V6.SOURCE instrument re- 
sponse functions without in-flight PSF corrections, 
selecting photon events between 70 MeV and 300 
GeV. Compared to the more typical 100 MeV 
threshold this choice leads to additional system- 
atic errors due to increased dispersion in the pho- 
ton energy reconstruction. The overall systematic 
flux error is energy dependent: it amounts to 30% 
at 70 MeV and dec reases to 10% above 10 GeV 
(jRando et al. 1 120111 ). The dominant part of this 



systematic error is related to the overall flux nor- 
malization. It is caused by uncertainties in the 
effective area determination and of the overall nor- 
malization of energy scale. 

The stability of the LAT instrument over time 
was tested for all time scales addressed in this 
publication using the Vela and Gem inga pulsars 



which are found to be stable in flux (Abdo et al 



2010bU cD. Their flux variations were <10%, yield- 



ing an upper limit on the variations of the sys- 
tematic errors with time. A more detailed study 
of the systematic uncertainties is currently being 
prepared within the LAT collaboration for publi- 
cation. 

The pulsar phase was assigned to the detected 
gamma rays based on a high-time-resolution pul- 
sar ephemeris. To obtain the latter we extracted 
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Fig. 1. — Phase profile of gamma rays above 70 
MeV within 3° of the Crab Pulsar for the first 33 
months of Fermi observations (black histogram) 
and during the April 2011 gamma-ray flare be- 
tween MJD 55663.70-55671.02 (black markers 
with error bars). The gray region indicates the 
adopted off-pulse interval where the emission is 
dominated by the nebula. The flare phase profile 
has been multiplied by a factor 59, such that the 
excess above the off-pulse counts is the same as for 
the 33 months of observations. This demonstrates 
that the flare is a phase-independent flux increase. 

400 pulsar times-of-arrival (TOAs) from LAT pho- 
tons collected from MJD 5 4684-55668 with a typi- 
cal uncertai nty of ^lOO^s (iRav et al. 2011 ). Us- 



ing tempo2 dHobbs et all ,20061) . we fit a timing 



solution to these TOAs with a typical residual of 
108^s, or about 3 x 10~ 3 of the pulsar period. To 
obtain these white residuals, we mo deled the pul 



sar tim ing noise using the method of lHobbs et al 
([20041 ) . with 20 harmonically related sinusoidal 
terms. The ephemeris parameter file, and light 
curves and spectra shown in this publication are 
publicly available onlin^E 

3. Time- Average Energy Spectra 

The Crab appears to the LAT as a point source, 
even at the highest photon energies. To sepa- 
rate nebular gamma-ray emission from that of the 
pulsar we apply a cut on the pulsar phase. The 
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Fig. 2. — Spectral energy distribution for the 
Crab nebula averaged over the first 33 months of 
Fermi observations. The axis on the right side 
indicates the isotropic luminosity. Also shown 
are data from COMPTEL in the soft gamma- 
ray band ( Kuiper et al. 200ll ) and very high en- 
ergy gam ma-ray measurements fr o m Cherenkov 
telescopes ( Aharonian et al. 20061 Albert et al 



20111 ). The dashed line shows the maximum like- 
lihood model in the parametrization described in 
the text. 

phased count rate of the pulsar with its double- 
peaked structure is shown in Figure [TJ The pul- 
sar dominates the phase-averaged gamma-ray flux, 
but its flux in th e off-pulse interval from 0.56 to 



0.88 is negligible (jAbdo et al. Il2010al ). It is in this 



interval that we measure the properties of the neb- 
ula. 

The LAT detects the nebula in the energy 
range between the the high-energy end of the syn- 
chrotron and the low-energy end of the inverse- 
Compton components of the SED. The average 
nebular spectrum measured during the first 33 
months of Fermi observations is shown in Fig- 
ure [2j We fitted it as the sum of synchrotron 
and inverse-Compton components. The differen- 
tial photon spectrum, <$>(E), of the synchrotron 
component was parametrized with a power law 



F s {is - 1) 
(100 MeV) 1 -^ 



E 



-7s 



(1) 



1 http: / /www-glast. stanford.edu/pub_data/691/ 



with an integral photon flux above 100 MeV of 
F s = (6.1±0.2) x 10~ 7 cm- 2 s" 1 and a soft photon 
index of 75 = 3.59 ± 0.07. 
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The spectrum of the inverse- Compton compo- 
nent softens significantly with respect to a power- 
law at higher energies, as expected from mea- 
sureme nts at very high energy (VHE) gamma 



rays ( Aharon ian et al. 1 120061: lAlbert et al. 1 12011 



lAbdo et al. Il2010al) . The spectrum was parametrized 



by a smoothly broken power law with a curvature 
index of /3 = 0.2: 




The spectrum has a photon index of 7^1 =1.48 
± 0.07 at low energies and softens to 77.2 = 
2.19 ± 0.17 above a break energy of E b = (13.9 
± 5.8) GeV. The flux normalization is $/,o = 
(5.6 ± 1.4) x lO" 10 cm~ 2 s" 1 MeV" 1 and the inte- 
gral flux above 100 MeV is F 7 = (1.1 ±0.1) x 10~ 7 

—2 —i 
cm s . 

The averaged pulsar spectrum in the first 33 
months of observations was measured in the on- 
pulse after accounting for the nebula emission. We 
parametrized the pulsar spectrum with a power- 
law function with a super-exponential cutoff 



E \ — 1P 
$p(E) = $p i0 I : - - I e 



1 GeV J 



(3) 



The best-fit value for the spectral index is 7p = 
1.59±0.01, for the break energy E P , C = (504±63) 
MeV and for the curvature index n — 0.43 ± 0.01. 
The normalization is given by $p.o = (8.1 ±0.5) x 
10~ 10 cm~ 2 s _1 MeV -1 . The integral flux above 
100 MeV is F P = (20.4 ± 0.1) x 10~ 7 cm' 2 s _1 , 
in agreement with the v alue measured after eight 



months of observations (jAbdo et al. 1120 10a') 



The pulsed emission from the Crab has recently 
been detected between «25-400 GeV. This emis- 
sion is constrained to pulses which are narrower 
by a factor of approximately two compared to the 
LAT energy range. Pulsed emission at these ener 



gies challenges current pulsar m odels (|Aliu et al. 



2008. 1201 it Uleksic et al~ll201ll) . A more detailed 



spectral analysis of the Crab pulsar in the Fermi- 
LAT energy range and its connection to the pulsed 
VHE emission is currently being performed and 
will be published elsewhere. 



4. Temporal flux variations 

In order to probe the flux variation over time, 
the flux from the Crab was evaluated in 12-hour 
time intervals. The combined pulsar and nebular 
spectrum was modeled as a power-law in energy, 
as further spectral features are not resolvable on 
these short time scales for typical observed fluxes. 

The resulting light curve is shown in Figure [3J 
where the three flares previously reported by the 
Fermi- LAT collaboration are indicated. Variabil- 
ity is seen over the whole 35-month period. This 
is also apparent in the Fourier power density spec- 
trum (PDS) shown in Figure 01 which is approx- 
imately described by a power law with a spectral 
index of «0.9. The PDS was derived using the fast 
Fourier transform for evenly spaced data, inter- 
polating over the short time interval MJD 54901- 
54906, during which the Crab was not observed by 
Fermi. We verified that results are essentially in- 
distinguishable from power spectra computed us- 
ing the techniques for unevenly sampled time se- 
ries described later in this section. The PDS shows 
significant power above the noise level on time 
scales from years to weeks. Variations on these 
time scales are also present outside of the flaring 
periods. This can be seen in the PDS of the time 
interval MJD 54884-55457 also shown in FigurelH 
where the nebula did not show large variations in 
flux. The spectral index of the PDS during this 
time is approximately 1.0. 

There are compelling reasons to believe that 
the phase-averaged Crab pulsar and the inverse- 
Compton component flux are constant i n the 



LAT energy range (|de Jager fc Harding I [1992; 



Atovan et al. 1 119961 ). or at least that they vary 
only very slowly as the pulsar spins down, and 
that any rapid variability is attributable to the 
synchrotron component of the nebula. This was 
confirmed observationally on monthly time scales 
for the first 24 mo nths of Fermi observations 
(|Abdo et al. Il2011ah . To test this assumption on 



the 12-hour time scales studied here, we measured 
the flux in the on-pulse interval, subtracting the 
nebula flux measured in the off-pulse interval for 
each time bin. The pulsar flux was found to be 
stable within 20%. No significant variations were 
found. The inverse-Compton component of the 
nebula is too faint to be consistently detected in 
12-hour time windows. Therefore the variability 
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Fig. 3. — Integral flux above 100 MeV from the Crab as a function of time for the first 35 months of Fermi 
observations. The time binning is 12 hours. The dotted blue line indicates the sum of the 33-month average 
fluxes from the inverse-Compton nebula and the pulsar. The dashed blue line shows the average flux of 
the synchrotron nebula summed to the latter. Fluxes are shown on a log scale. The three flares detected 
by the Fermi L AT are indicated. A fourth flare was detected in October 2007 (MJD 54380) by AGILE 
(jTavani et al. 11201 ll ). 
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Fig. 4. — Crab nebula Fourier power density spec- 
trum (PDS), calculated from the light curve of 
the first 35 months of Fermi observations shown 
in Figure [3J The PDS of the full time interval 
is shown by the solid green line (scaled down by 
1/100 for better visibility). The PDS of the low ac- 
tivity period between MJD 54884-55457 is shown 
by the solid blue line. The PDS of the April 2011 
flare is indicated by the solid red line and was 
calculated from the light curve shown in Figure 
[5l A smoothing with a running average of four 
bins was applied to all spectra. The PDSs ob- 
tained before smoothing are shown in colored dot- 
ted lines. Black lines show the best fit function of a 
power-law function (dashed) plus a constant white 
noise component (solid) for the unsmoothed spec- 
tra. The best fit spectral indices are given in the 
text. Dotted black lines indicate the ±lcr, the +2<t 
and +3<7 confidence intervals derived from white 
noise simulations for the April 2011 flare PDS. 



of this component cannot be quantified from our 
analysis on these time scales. However, significant 
flux variations of the inverse- Compton component 
on short time scales are strongly disfavored the- 
oretically. The observed absence of variations on 
monthly time scales is consistent with constancy 
on the shorter scales of interest here. 

On 2011 April 9 the flux of the Crab increased 
dramatically. Once the source reached flux levels 



comparable to the 2009 and 2010 flares the Fermi 
satellite was commanded to switch to a pointed 
observation targeting the Crab (MJD 55663.70- 
55671.02). In this mode the exposure toward the 
source increased by a factor of about four com- 
pared to the standard all-sky monitoring. During 
these observations the Crab erupted to a peak flux 
of (186 ± 6) x 10~ 7 cm" 2 s^ 1 above 100 MeV dur- 
ing the 12-hour period centered at MJD 55667.14, 
as shown Figure [3J This corresponds to a flux 
increase of a factor 7 compared to the average to- 
tal flux from the Crab and a factor of 30 com- 
pared to the flux from the nebula. We verified that 
the flare is indeed positionally coincident with the 
Crab nebula, with a best-fit localization of R.A. = 
83.65° Dec. = 21.98° (J2000), and an error radius 
of 0.04° at 95% confidence. 

The high gamma-ray flux during the flare, com- 
bined with the increased exposure of the pointed 
observations allowed us to study the flux evolution 
down to time scales below the wl.5 hour orbital 
period of Fermi. On such time scales occultation 
of the source by the Earth needs to be taken into 
account. The time binning was therefore adjusted 
such that only time periods during which the Crab 
is visible to the LAT are used. These visibility win- 
dows were further split into bins of equal exposure, 
yielding a mean bin duration of nine minutes. The 
evolution of the flux during the flare in this binning 
is shown in Figure [5] The flare lasted for approx- 
imately 9 days and is composed of two sub-flares, 
peaking around MJD 55665 and MJD 55667. Dur- 
ing both, the flux increased rapidly, reaching its 
maximum value within approximately one day. A 
second rise or "shoulder" is observed during the 
decaying phase of both sub-flares. Whether this 
is coincidental cannot be assessed on the basis of 
these two events alone. If these shoulders are not 
interpreted as additional flares of lower amplitude, 
the decay time is approximately 1.5 and 3 days for 
the first and second sub-flare, respectively. 

To investigate the significance of peaks on 
smaller time scales in the flare light curve, we 
decomposed it into periods compatible with a con- 
stant flux. The time windows are referred to as 
"Bayesian Blocks" (BBs). They were determined 
finding the partition which maximizes th e sum of 



the c ost function assigned to each BB (jScargle 



19981 ). The cost function was set to the loga- 



rithm of the maximum likelihood for the constant 
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Fig. 5. — Integral flux above 100 MeV as a function of time during the 2011 April Crab flare. The light 
curve is binned into equal exposure bins during times with no Earth occultation, with a mean bin duration 
of nine minutes. The dotted line indicates the sum of the 33-month average fluxes from the inverse- Compton 
nebula and the pulsar. The dashed line shows the flux of the average synchrotron nebula summed to the 
latter. The solid black lines show the best fit of a model consisting of a constant plus an exponential function 
at the rise of both sub-flares (see text). The blue vertical lines indicate the intervals of each Bayesian Block 
during which the flux remains constant within statistical uncertainties. The time windows are enumerated 
at the top of the panel. The corresponding flux is shown by the blue marker below each number. The SED 
for each of the time windows is shown in Figure [6] 



flux hypothesis. The algorithm to determine the 
optim al partition is described by Ijackson et al. 
(|2005h . The BB-binned light curve is shown in 
Figure [5] It is statistically compatible with the 
original light curve {x^/ndf = 257/232). This im- 
plies that flux variations within each BB cannot be 
distinguished with confidence from a locally con- 
stant flux. The shortest BBs are detected at the 
maximum of both sub-flares and have durations 
of «9 hours. 

In order to measure the rate of flux increase at 
the rising edges of the sub-flares we parametrized 
them with an exponential function plus a constant 
background. The best-fit functions are shown in 
Figure [5] The time ranges over which the fits were 
performed were defined by the centers of the BBs 
before and at the maximum of each sub-flare. The 
resulting doubling time is 4.0 ±1.0 hours and 7.0 
±1.6 hours for the first and second sub-flare, re- 
spectively. As these values depend on the some- 
what arbitrarily chosen parametrization and fit 



ranges, we conservatively estimate that the dou- 
bling time scale in both sub- flares is t<j < 8 hours. 

The PDS of the April 2011 flare is shown in 
Figure [4] It was obtained by computing the 
Fourier transform of the autocorrelation function 
using an algorithm for unevenly sampled data 



(|Edelson fc Krolik I Il988h . The PDS can be de- 
scribed by a power law of index wl.l and reaches 
the noise floor at a frequency of «0.6 cycles per 
day. The doubling time of the corresponding sinu- 
soidal component is ~10 hours, in agreement with 
the expectation from the measured doubling times 
of the flares. 

The pulsar flux remained unchanged during 
the flare, with an average flux above 100 MeV 
of F P = (21.7 ± 1.1) x 10~ 7 cm- 2 s" 1 dur- 
ing the main part of the flare (MJD 55663.70- 
55671.02). The flux increase is phase-independent. 
This is illustrated in Figure [TJ where the pha- 
sogram during the main flare period is shown. 
The peaks in the on-pulse interval remain at the 
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Fig. 6. — Spectral Energy Distribution evolution during the April 2011 Crab flare. Arrows indicate 95% 
confidence flux upper limits. The time windows are indicated in the bottom left corner of each panel and 
correspond to the ones indicated in Figure [5j The dotted line shows the SED of the flaring component, the 
dot-dashed line the constant background from the synchrotron nebula, and the dashed line is the sum of both 
components (see text). The average Crab nebular spectrum in the first 33 months of Fermi observations is 
also shown in gray for comparison. 



same position. We also searched for periodici- 
ties other than the Cr ab pulsar with the time- 



differencing technique (jAtwood et al. 2006) 



ap- 



plying the ev e nt-we ighting technique described in 



Bickel et al. I ([20081 ) . We scanned the frequency 
range 0.1-256 Hz, allowing for a possible spindown 
up to twice the value of the Crab pulsar. No sig- 
nificant signal was found besides the pulsar, which 
was detected with a significance > 5.5er. Finally, 
we searched for photon clumping on time scales 
shorter than the «10 min time binning by apply- 
ing a Bayesian Block analysis on the single photon 
arrival times, with no significant detection. 



4.1. Spectral evolution during the flare 

In order to measure the energy spectrum dur- 
ing the flare, and its evolution with time, the data 
must be averaged in time intervals long enough 
to ensure adequate photon statistics, but short 
enough to provide adequate temporal resolution. 
The 11 bins of approximately constant flux, de- 
rived from the BB analysis, provide a reasonable 
compromise between these two constraints. 

The SEDs for each of the time bins are shown 
in Figure El after subtracting the steady emis- 
sion from the pulsar and the inverse- Compton 
component of the nebula. It can be clearly seen 
that a new spectral component emerges from the 
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synchrotron nebula during the flare, moving into 
the Fermi energy range as the flare evolves. Its 
flux reaches a maximum between MJD 55666.997- 
55667.366 (frame 7); during this period the peak 
in the SED is clearly detected at E pea k = (375 ± 
26) MeV. 

It is difficult to parametrize the spectral shape 
of the flaring component due to the likely contam- 
ination by background flux from the synchrotron 
nebula not related to the flare. The determination 
of the latter is degenerate with the measurement 
of the flare component, as only the summed flux is 
measured. To break this degeneracy we proceeded 
under the following assumptions: 

1. The spectrum of the synchrotron nebula 
during the flare can be described by a power- 
law function and does not vary in time. 

2. The spectrum of the flaring component can 
be described by a power-law with an ex- 
ponential cutoff (equation [3] with k, = I). 
While the cutoff energy and normalization 
of the spectrum vary, the spectral index re- 
mains constant during the flare. 

We derived the spectral index of the flaring com- 
ponent and the spectrum of the background syn- 
chrotron component in a composite likelihood fit 
to all the time windows displayed in Figure [5j si- 
multaneously measuring the energy cutoff and flux 
normalization evolution of the flaring component 
in each of the time windows. For this we used 
the composite likelihood 2 part of the Fermi 
Science Tools. 

The best-fit values for the background syn- 
chrotron nebula during the flare period are F$ = 
(5.4 ± 5.2) lCT 7 cm" 2 s" 1 and 7s = 3.9 ± 1.3, 
consistent with the average value measured during 
the first 33 months of observation. The spectral 
index of the flaring component is measured to be 
7_f = 1.27 ± 0.12. The best-fit values for Ep tC and 
the energy flux above 100 MeV, fp, are shown in 
Figure [7J 

This simple parametrization yields a good de- 
scription of the flare evolution, as can be seen in 
Figure |H1 The data are therefore compatible with 
the interpretation that a new spectral component 
of a fixed spectral shape evolves during the flare. 
As shown in Figure [7] the cutoff energy of the flar- 
ing component varies significantly, having a statis- 
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Fig. 7. — Total energy flux above 100 MeV, fp, as 
a function of cutoff energy E c for the flaring com- 
ponent of the April 2011 flare. The values corre- 
spond to the parametrizations shown in frames 1— 

10 of Figure [6] (the values obtained from frame Nr 

11 are not included, as no significant spectral cur- 
vature was detected in his time interval, allowing 
no robust determination of E c ). The number next 
to each marker denotes the corresponding frame. 
The axis on the right hand side indicates the ap- 
parent luminosity in units of the pulsar spin -down 



power of 5 x 10 38 ergs s _1 (|Hester I [200^) . The 



numbered solid lines indicated the slope of the 
corresponding power-law dependency fp ~ Ep c a . 
The dotted line marks the best fit function in this 
parametrization with a = 3.42 ± 0.86. 

tical probability of being constant of only 0.04%. 
However, while fp varies by more than an order 
of magnitude, the cutoff energy varies only by ap- 
proximately a factor of two. The relationship be- 
tween both quantities can be described approxi- 
mately by a power-law function fp ~ Ep^ c a , with 
best-fit index of a = 3.42 ± 0.86. 

5. Discussion 

A year after first being reported, the gamma- 
ray flares from the Crab remain enigmatic. Where 
within the nebula does the emission come from? 
What produces the flux variations? How were the 
emitting particles accelerated? How are the flares 
related to the variability observed on yearly and 
monthly time scales? Although several ideas have 
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been proposed, no certain answers can be given 
today. The observations presented here give us 
the most precise look into the flare phenomenon to 
date, during the brightest outburst detected so far. 
We will proceed to discuss some of the implications 
and challenges posed by these observations. 

One striking property of the Crab nebula flares 
is their rapid flux variations, doubling within t^ < 
8 hours at the rise of the 2011 April flare. Causal- 
ity arguments imply that the emission region is 
compact, with a length L < ctd ~ 2.8 x 10~ 4 
pc. The emitted isotropic power at the peak of 
the flare of « 4 x 10 36 erg s _1 corresponds to 
wl% of the total spin-down power of the pul- 
sar, the ultimate energy source of the nebula. 
It is difficult to explain how this energy is fo- 
cused into such a small emission volume. The 
focusing is generally easier to explain when the 
emission site is closer to the pulsar. The ab- 
sence of pulsation in the flare signal implies that 
the emission region is at least located outside the 
light cylinder of the pulsar. Another possibility 
to explain the flare brightness is that the emis- 
sion is highly anisotropic, as would be expected if 
the emission region moves relativistically toward 
us. While only mildly relativistic motion with ve- 
l ocities of ~0 . 5c are observed inside the nebul a 
( Scarglelll969llHester et al. Il2002t [Melatos et al 



relativistic motion is expected in the pul 
sar wind and in the downstream medium behind 



the wind termination ( Camus et al. 1 l2009h . Rel- 
ativistic bulk motion is particularly expected at 
the "arch shock" of the wind termination, which 
has been proposed as the main site of ga mma-ray 
emission (|Komissarov fc Lvutikov 112011 ). 

The flare emission is expected to result from 
synchrotron radiation by relativistic electrons and 
positro ns (henceforth refer red to together as elec- 



trons) (jAbdo et al. 11201 laf ). A new spectral com- 



ponent emerges in the SED during the flare. The 
hard photon spectrum (7 » 1.3) of the flaring 
component implies that most of the electron en- 
ergy is carried by the highest energy electrons. 
If the electron particle density n(e) per energy 
at an energy e is characterized by a power law 
n(e) ~ e~ p , the spectral index is p < 2-f — 1 » 1.6, 
in a random magnetic field in which the elec- 
trons are isotropically distributed. The energy 
per logarithmic energy interval, which is propor- 
tional to e 2 x ?i(e), is therefore rising with en- 



ergy. We note that such a spectrum is also in- 
ferred for the radio emitting electrons in the Crab 
neb ula, and more generally in pulsar wind nebu- 
las (Gaensler fc Slane 2006 ; Sironi &: Spitkovskv I 



2011al ). While the radio emission is produced by a 
different electron population than the gamma-ray 
flares, efficient particle acceleration appears to be 
a common feature in these systems. 

It is an interesting question how such a hard 
electron spectrum is produced. Standard diffu- 
sive shock accel eration typically results in spec 



tra w ith p > 2 ([Gallant et al. 1119921 iKirk et al. 
2000). Even though it has been shown that harder 



spectra can be produced in certa in field configura- 
tions with low-level turbulence (|Kirk fc Heavens! 



1989; Summerlin fc Baring 2011 ). these condi 



tions are not expected at termination shocks of 
pulsar winds. Additionally, shock acceleration ap- 
pears to be inefficient at highly oblique shocks 
that are representative of the pulsar wind ter- 
mination discontinuity (Ellison fc Double 2004; 



Summerl in &: Baring 1 1201 It ISironi &: Spitkovskv 



2011bl ). One alternative is that magnetic re- 



celerate particles ( 


Lvubarskv 2003 


: Kirk 


2004; 


Yuanetal. 2011; 


Bednarek et al. 1 


2011). 


How- 



ever, simulations show that reconnection behind 
the pulsar termination probably does not provide 
the required e lectron energies to produce g amma- 
ray emission ([Sironi fc Spitkovskv I l2011al ) . An- 
other interesting possibility is that acceleration is 
occurring directly in the e lectric field induced b y 
the pulsar, as discussed by lAbdo et al. 1 (|2011al ). 

The observation of a peak synchrotron energy 
of ~380 MeV is among the highest yet seen from 
astrophysical source today. The observation is 
surprising as particle acceleration in the presence 
of synchrotron cooling is expected to limit syn- 
chrotr on emission to photon energies below ~ 150 
MeV dGuilbert et al. Ill983t Ide Jager et al. I fl996l: 



Komissarov fc Lvutikov 20 111) . Two solutions to 
this problem have been proposed recently in this 
context: 

1. The electric field at the acceleration site 
is larger in magnitude than the magnetic 
field. This is generally an unstable state 
in plasma, as charges will short out the 
electric field; however, temporarily such 
a configuration is expected, e.g. in mag- 
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netic reconnection events. The gamma- 
ray emission might occur afterward, when 
the accelerated electrons enter a region of 
enhanced magnetic fields (lUzdenskv et al. 



201ltlCerutti fe Uzdenskv Il201l[) . More de- 



tailed studies are required to assess whether 
such a scenario is plausible in the nebula 
environment and can be sustained for the 
duration of the flares. 

2. The gamma rays are emitted in a region 
of bulk relativistic motion, and are there- 
fore Doppler boosted toward the observer 
(jKomissarov fe Lvutikov 201 lh . For a flow 
moving directly toward us a Lorentz factor 
> 2 is sufficient to accommodate the ob- 
served peak energy. 

Variations in the Doppler boosting can nat- 
urally account for th e observed flux variation 



(jLvutikov et al. 1 1201 lh . The observed spectral 



evolution is compatible with such an interpreta- 
tion: the energy flux of the emission varies ap- 
proximately as a power of a = 3.42 ± 0.86 with 
the cutoff energy; a correlation with a » 3 is in- 
deed expected for varia tions produced by change s 
in relativistic beaming ([Lind fe Blandford 1119851 ). 

The flare brightness, the high frequency of the 
observed peak of the gamma-ray emission, and the 
spectral evolution during the flare all suggest the 
presence of relativistic beaming. We therefore con- 
clude that, independent of the location of the emis- 
sion region and the physical processes responsible 
for the flares, the emission region is moving rela- 
tivistically toward us, and changes in its motion 
are likely the predominant mechanism responsi- 
ble for the observed flux variations. Such a kine- 
matic explanation does however not address the 
issue of how a moving source can be created dy- 
namically and sustained radiatively in the face of 
strong losses. The Crab Nebula still has much 
more to teach us. 
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